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We have found that at least seven hydrogen-deficient carbon (HdC) and R 
Coronae Borealis (RCB) stars, have 16 0/ 18 ratios close to and in some cases 
less than unity, values that are orders of magnitude lower than measured in other 
stars (the Solar value is 500). Greatly enhanced 18 is evident in every HdC and 
RCB we have measured that is cool enough to have detectable CO bands. The 
three HdC stars measured have 16 0/ 18 < 1, lower values than any of the RCB 
iy-j \ stars. These discoveries are important clues in determining the evolutionary 

pathways of HdC and RCB stars, for which two models have been proposed: 
the double degenerate (white dwarf (WD) merger), and the final helium-shell 
flash (FF). No overproduction of 18 is expected in the FF scenario. We have 
quantitatively explored the idea that HdC and RCB stars originate in the mergers 
of CO- and He-WDs. The merger process is estimated to take only a few days, 
with accretion rates of 150 M yr~ l producing temperatures at the base of the 
accreted envelope of 1.2 — 1.9xl0 8 K. Analysis of a simplified one-zone calculation 
shows that nucleosynthesis in the dynamically accreting material may provide a 
suitable environment for a significant production of 18 0, leading to very low values 
of 16 0/ 18 0, similar to those observed. We also find qualitative agreement with 
observed values of 12 C/ 13 C and with the CNO elemental ratios. H-admixture 
during the accretion process from the small H-rich C/O WD envelope may play 
an important role in producing the observed abundances. Overall our analysis 
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shows that WD mergers may very well be the progenitors of 18 -rich RCB and 
HdC stars, and that more detailed simulations and modeling are justified. 

Subject headings: stars: evolution — stars: abundances — stars: AGB and post- 
AGB — stars: variables: R Coronae Borealis 



Introduction 



Among the hydrogen-deficient post-asymptotic giant branch (post- AGB) stars are the R 
Coronae Borealis (RCB) stars, a small group of carbon -rich supergiants. About 50 RCB stars 



are known in the Galaxy and the Mage llanic Clouds (IKraemer et al.l 120051 ; IZaniewski et al. 



20051 : lAlcock et al.l l200ll ; IClaytonl Il996l ). Their defining characteristics are hydrogen defi- 
ciency and unusual variability - RCB stars undergo massive declines of up to 8 mag due to 
the formation of carbon dust at irregular intervals. Apparently related to the RCB stars 
are the hydrogen-deficient carbon (HdC) stars. The five known HdC s tars are simil ar to the 
RCB stars spectroscopically but do not show declines or IR excesses (jWarnerlll967l ). 



Two scenarios have been proposed for the origin of an RCB star: the double degenerat e 
(DD) and the final helium-shell flash (FF) models jlben et ali[l996: ISaio fc Jefferv! 120021 ) . 
The former involves the merger of a CO- and a He-white dwarf fWD) JWebbmkHl984h . In the 
latter, thought to occur in 20% of all AGB stars, a sta r evolving into a planetary nebu la (PN) 
central star is blown up to supergiant size by a FF (jFujimotol 119771 ; lRenzinilll979l ). Three 
stars (Sakurai's Object, V605 Aql, and FG Sge) have been observed to undergo FF outbursts 
that transformed them from h ot evolved PN central stars into cool giants with spectral 
properties similar to RCB stars JAsplund et alill998l . Il999l . I2OO0I : Idavton fc De Marcdfl997l : 
Gonzalez et al.|[l998 ). In two of these, Sakurai's Object and V605 Aql, the FFs are thought to 



be th e result of Ver y Late Thermal Pulses (VLTPsjl, FFs that occur on the WD cooling track 



Herwig!l200ll ). because their post-outburst evolutions have been so rapid and be cause 



(e.g. 

in some cases an ionized nebula was found surrounding the cool post-outburst object (IFord 
197ll : buerbeck fc Benett] Il996h . V605 Aql has evolved from T e// =5000 K in 1921 (two 
years after its outburst) to 95,000 K today Jdavton et alil2006h and now has abundances 
similar to those seen in Wolf-Rayet [WC] central stars of PNe, with 55% He, and 40% C, and 
not similar to those seen in the RCB and HdC stars. Sakurai's Object became enshrouded in 
dust three years after its outburst, and thus most of its current stellar properties cannot be 



1 FG Sge is thought to be the result of a Late Thermal Pulse (LTP) and is not considered further here 
(|HerwiEll200lT ). 
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determined. However, just after their outbursts, both V605 Aql and Sakurai's Object were 
almost indistinguishable from the RCB stars, in abundances (except for the presence of 13 C), 
temperature, and absolute luminosity. Nevertheless, the very short observed timescales for 
this RCB-like phase make it unlikely that object s such as these can account for even the 
small number of RCB stars known in the Galaxy (jClayton et al.ll2006l ). 



Recently, IClavton et all (120051 ) discovered that 16 0/ 18 < 1 in the HdC star, HD 
137613, from a spectrum of the first overtone bands of CO at 2.3-2.4 /zm. As the bands of 
i2q16q are Q f typi ca i strength for a cool star, the measured ratio implies a huge enhancement 
of 18 O rather than a depletion of 16 O. The isotop ic ratio is orders of magnitude less than 
that measured in ot her stars (|Clavtqn et al.ll2005l). Th e value of 16 0/ 18 is ~500 in the 



solar neighborhood (jScott et al 



2006 



Gciss et al. 



20021 ) and varies from 200 to 600 in the 



Galactic interstellar medium ( Wilson fc Roodl " 1994 ). Only one other star, the post- AGB 



star HR 4049, has been found to have highly enhanced 18 (ICami fc Yamamurall200ll ). HR 
4049 is a binary and its enhanced oxygen isotopes ( 16 0/ 18 ~ 7 and 16 0/ 17 ~ 8) are 
found in the circumbinary disk material. High abunda nces of 18 O have al so been measured in 
pre-solar graphite grains in the Murchison meteorite (lAmari et al.lll995l ) and are attributed 
to material processed in massive Wolf-Rayet stars. In both of these other cases (HR 4049 and 
meteorites), the abundance of 18 relative to 16 is several times less than in HD 137613. 

In this paper, we analyze new high S/N spectra of the K-band CO bands in HD 137613 
and nine additional RCB and HdC stars to investigate 16 0/ 18 in these stars. The sample 
includes all known H dC stars and five RCB stars cool enough to have prominent CO bands. 
Clayton et al.l (120051 ) discussed the failure of several nucleosynthesis scenarios to explain the 
unusually small 16 0/ 18 ratio in HD 137613. In the light of the new observations, we explore 
the idea that HdC and RCB stars originate in DD mergers of a CO and a He- WD. We first 
investigate whether the conditions exist in a WD merger for nucleosynthesis to take place 
and then estimate the results of the nucleosynthesis. 



2. Observations and Data Reduction 



We obtained a more extensive K-band spectrum of HD 137613 than the one published 
by IClavton et all fj2005h on UT 2005 March 04 at the 3.8 m United Kingdom Infrared 
Telescope (UKIRT). The new observation used the facility instrument UIST with its long 
K grism and Of! 24 (2 pixel- wide) slit to cover 2.20-2.51 fim at a resolving power of ~3500 
( IRamsay Howat et al.l l2004f ) . Standard observing and reduction techniques were employed. 
A wavelength calibration accuracy of better than 0.0001 fim was achieved using the spectrum 
of an argon arc lamp. The calibration star used to remove telluric absorption lines and flux- 
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calibrate the spectrum of HD 137613 was HR 5514 (B9V). The reduced spectrum is plotted 
in Fig. 1, with vertical lines indicating the bandhead wavelengths. I t show s, in addition 
to the 2-0, 3-1, and 4-2 bands of 12 C 18 discovered bv Idavton et all (b()05h . the 5-3, 6-4, 
and 7-5 bands. In addition to the six bands of 12 C 18 0, at least seven and probably eight 
bands of 12 C 16 are seen. As in the previously published K-band spectrum of HD 137613 
the strengths of the 12 C 18 and 12 C 16 bands are approximately equal, implying near equal 
abundances of the two isotopomers. The spectrum shows no evidence for either 13 C 16 or 

12 C 17 . 

The Gemini-S facility spectrometer, GNIRS was employed on UT 2005 September 26- 
28 to obtain spectra of nine HdC and RCB stars. The spectrograph was used in the short 
camera mode and with the 0'.'3 (2 pixel-wide) slit to cover 2.27 to 2.51 fim at R=5900. The 
observations were reduced by dividing the spectrum of each RCB and HdC star by that of a 
nearby bright early A-type dwarf star to remove the instrumental response and the effects of 
atmospheric absorption lines. Wavelength calibration was obtained from these telluric lines 
and is accurate to better than 0.0001 /jm. 

The complete observing sample is listed in Table 1 and the 2.28-2.45 /zm portions of 
the spectra are plot ted in Fig. 2. Also inclu ded is the previously published spectrum of the 
RCB star, Z UMi f lTenenbaum et al.l 120051 ) . The wavelengths of the bandheads of 12 C 16 
and 12 C 18 are indicated by vertical lines. The signal-to-noise ratios of the spectra are 100 
or more. The spectra have been shifted to correct for the stars' measured he liocentric radial 
velocities (ILawson Cottrelll 119971 ; iDuflot et al.lll995l ; iGoswami et al.lll997l ). One star, ES 
Aql, does not have a published radial velocity. We have estimated it to be ~55 km s _1 and 
made the correction in Fig. 2. 



3. Preliminary Results 

Seven of the eleven stars in Fig. 2 show, in addition to prominent bands of 12 C 16 0, com- 
pelling evidence for bands of 12 C 18 0. These are the HdC stars HD 137613, HD 175893, and 
HD 182040, and the RCB stars WX CrA, ES Aql, S Aps, and SV Sge. One additional RCB 
star that we observed, U Aqr, shows marginal evidence for 12 C 18 Q. We also find margin al 



evidence for 12 C 18 in the the spectrum of the RCB star, Z UMi (ITenenbaum et al.ll2005l ). 



In the spectra of three of the HdC stars, HD 137613 and HD 175893, and HD 182040, and 
in the spectrum of one RCB star, WX CrA, the 12 C 18 bands are about as strong or stronger 
than the 12 C 16 bands. The most remarkable case is the HdC star HD 175893, for which 
the maximum depths of the bands of 12 C ls O are more than twice as deep as those of 12 C 16 0. 
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The least secure of these four cases is the HdC star HD 182040, because the bands are weak, 
but there is accurate wavelength correspondence between observed absorption maxima and 
their expected wavelengths. As in the case of HD 175893, the 12 C 18 bands in its spectrum 
are noticeably stronger than the 12 C 16 bands. For the third HdC star, HD 137613, the 
i2q18q "bands appear to be slightly stronger than the 12 C 16 bands, whereas in the RCB 
star with the strongest 12 C 18 0, WX CrA, the band strengths do not exceed those of 12 C 16 0. 
In three other RCB stars, ES Aql, S Aps, and SV Sge, the 12 C 18 bands are prominent, 
but are considerably weaker than than the 12 C 16 bands. Two HdC stars, HD 148839 
and HD 173409, have no significant CO bands, presumably because their relatively high 
photospheric temperatures preclude the formation of detectable amounts of CO. Thus for 
them no determinations of the relative abundances of the CO isotopomers are possible. 

For the sources with weaker absorptions (U Aqr, Z UMi, and HD 182040), there is 
some possibility that the features at the expected wavelengths of the 12 C 18 bandheads 
are due to some other species (e.g., CN). However, the accurate wavelength agreement of 
the 12 C 18 bandheads (in terms both of absolute wavelength and wavelength difference 
from nearby 12 C 16 bandheads), plus the fact that the absorptions at those wavelengths 
are among the strongest (apart from the 12 C 16 absorptions) in each spectrum gives us 
confidence, especially in the case of HD 182040, that we are detecting 12 C 18 0. None of the 
spectra show evidence for absorption bands of 13 C 16 0, which are present in the spectra of 
some carbon stars. Very high 12 C/ 13 C ratios are one o f the trademarks of the RCB and HdC 
stars jAlcock et alihoOll ; lciavt~onlll996l : Iwarnerlll967h . 



4. Estimates of 16 0/ 18 



To further confirm the identification of the 12 C 18 bands in our sample and to provide 
a way to roughly estimate 16 Q/ 18 0, we have computed m arcs model atmosphere spectra for 
different CO isotopic abundances (lAsplund et al.lll997l ). We ran a grid of models with the 
stellar parameters T e ff = 5000, 5500, 6000 K and log g = 1.5 [cgs] and a chemical composition 
typical of RCB stars JAsplund et al.lll997h . Values of 16 0/ 18 chosen were, 0, 0.1, 0.33, 1, 
3, 10, and infinity (no 18 0). Models with no CO were also run. 

The spectra in Figs. 1 and 2 and the marcs model spectra are complex, because, 
in addition to CO, other species that are ab undant in carbon sta rs absorb strongly in the 
observed wavelength interval, as discussed by lClayton et al.l (120051 ). In order to estimate the 
16 0/ 18 ratio for each star we first need to measure the absorption depths of the 12 C 16 
and 12 C 18 2-0 and 3-1 bands. To do this, we started by binning the spectra in Fig. 2 and 
the marcs spectra in intervals of 0.0005 /im. From the marcs spectra with no CO, we then 
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identified wavelength bins in the "continuum" just shortward and longward of each bandhead 
where absorption by other species is relatively weak. The ratios of the flux densities at the 
bandheads to the average of those at the "continuum" wavelengths, corrected for the net 
absorptions at the bandheads in the marcs spectrum with no CO, are shown in the fourth 
column of Table 2 and provide a first estimate of the isotopic ratio 16 0/ 18 0. The values 
shown are probably overestimates of this ratio where 16 0/ 18 < 1 and underestimates where 
16 0/ 18 > 1, because the bands of the more abundant isotopomer are more saturated. 

We have interpolated between the marcs spectra to estimate the 16 0/ 18 ratio for each 
star. The values are given in the right hand column of Table 2, and confirm the very low 
ratios in seven stars. However, the complexities of the observed and model spectra make 
determination of the isotopic ratios problematical, especially in cases such as this, where 
individual lines are not resolved. The marcs spectra do not match the observed spectra 
in detail, adding to the uncertainties. We believe that the derived strengths of the 2-0 
bandhead of 12 C 18 are more reliable than those of the 3-1 bandhead, where the correction 
for contaminating species was largest, and so we have weighted the ratios of the 2-0 band 
strengths more highly. Finally, the depths of the CO features in S Aps and SV Sge are 
greater than those in any of the model spectra that we generated; for these two stars, we 
used the 5000 K marcs models. For all of these reasons we estimate the uncertainties in 
the derived isotopic ratios to be accurate only to a factor of two. 



5. Discussion 



In a prescient paper, IWarnerl (119671 ) predicted that HdC stars would be rich in 18 O. 
He noted that just before He burning begins in the core of an intermediate mass star the 
temperature increases to values where u N(a, 7) 18 F(/5 + z/) 18 can occur and subsequently the 



outer envelope of the star is lost throug h nova- or PN-like mass-loss episodes, exposing the 
18 0-rich material. Clayton et al. (|2005l ). armed with the detection of a high 18 O abundance 



in one HdC star, HD 137613, and using current information about nuclear reaction rates 
and stellar interiors, discussed the nucleosynthesis that could possibly lead to the unusual 
oxygen and carbon isotopic abundance ratios in an isolated intermediate mass star. They 
noted that the only evolutionary phase that is able to produce the observed oxygen isotopic 
ratio as well as satisfy the elemental abundance constraints for He and CNO, occurs at the 
onset of He-burning after completion of H-burning, via the same a-capture reaction as above. 
For example, in e arly-AGB star s the 18 Q enhancement is found only in a thin layer of the 
He-burning shell. IClayton et al.l (120051 ) speculated that HD 137613 may be a peculiar object 
that has somehow managed to strip all of its envelope material precisely down to the narrow 
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layer where the observed extremely low value of O/ O is found 



However, the observations reported here demonstrate that an extreme excess of 18 
is a common and possibly universal feature of HdC and RCB stars. Even if the above a- 
capture reaction is the source of the enhanced 18 O, the idea that for a large set of stars 
post-AGB mass loss proceeds nearly precisely to the 18 0-enriched layer and then stops is 
highly improbable. Another mechanism must be the cause of the enhancement. Moreover, 
the uniquely low values of 16 0/ 18 in HdC and RCB stars make it obvious that these two 
classes of carbon-rich and hydrogen-poor objects are indeed closely related. Therefore, we 
suspect that the isotopic ratios are an important clue in determining the formation pathway 
of HdC and RCB stars in general. 

Recent studies of V605 Aql an d Sakurai's Object i mply that a FF origin for most of 
the RCB and HdC stars is unlikely (jClayton et al.l 120061 ). based, for example, on the short 
time that FF stars spend in the region of the HR diagram where RCB stars are found. 
Although Sakurai's Object, as a certain representative of the FF scenario, s hares a number 



of characteristics with the RCB stars, there are some important differences (jAsplund et al. 



998. 



19991. 200 0j). Its abundances resemble V854 Cen and the other minority-class RCB stars 



(jAsplund et al 



19981 ). However, in Sakurai's Object and probably in V605 Aql, C/ C 



is low indicating efficient H-burning, while in the RCB and HdC stars it is generally high, 
thereby excluding th e possibility that we are seeing material processed by efficient H-burning 
(jClayton et al.ll2006l ). Moreover, infrared spectra of Sakurai's Object obtaine d during 1997- 
1998 when it had strong CO overtone bands showed no evidence for 12 C 18 (IGeballe et al. 



20021 ). Therefore, Sakurai's Object and the other FF stars on the one hand, and most of the 
RCB and HdC stars on the other hand, are likely to be stars with different origins. 

As we have discussed before, 18 can be overproduced in an environment of partial 
He-burning in which the temperature and the duration of nucleosynthesis are such that the 
1A N(a, 7) 18 F(/? + ) 18 reaction chain can produc e 18 Q, if the 18 Q is not further processed by 
18 0(a, 7) 22 Ne. This has also been suggested by Lambert (1986)0. Current understanding of 
mixing and nucleosynthesis events in the FF scenario, which leads to predictions in excellent 
agreement with Sakurai's Object (lAlcock et al.ll200ll ). is that the surface abundance after 
the FF is dominated by the intershell material, possibly with the superposed signature of H- 
ingestion nucleosynthesis if the FF was a VLTP. In the He-shell flash 14 N is completely burned 
into 22 Ne as temperatures are so high that 18 is destroyed. In the case of a VLTP, however, 



2 However, it should be noted that two of th e RCB stars, V3795 Sgr and Y Mus, as well as the extreme 
helium (EHe) stars, show high [Ne/Fe] ratios ( Asplund et al. 2000( l . These stars are too hot to show CO 
bands so nothing is known about their 18 abundances. 
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there is additional nucleosynthesis due to the ingestion of H-rich envelope material into the 
He-shell flash convection zone. This may separate the upper H-flash driven convection zone 
from the original He-shell flash convection zone. In the upper convection zone, He-burning 
will be terminated and one needs to consider if conditions for 18 overproduction exist 
there. However, p-capture nucleosynthesis, induced by the H-ingestion event, efficiently 
destroys 18 O and therefore we do not predict any overabundance of 18 O in the FF case. 
It should be mentioned though that the FF scenario involves complicated, hydrodynamic 
mixing events that may violate some of the assumptions made in 1-D stellar evolution. We 
cannot, therefore, conclusively settle this question at present. 

The above difficulties in producing an extreme 18 excess together with concerns about 
the ability of FF events to account for the observed number of RCB and HdC stars have led 
us to explore an alternative evolutionary origin for the RCB and HdC stars, namely that of 
DD mergers. In particular we have c o nsider ed the onset of He-burning nucleosynthesis in 
a He- WD/CO- WD merger. IWebbinkl (119841 ) proposed that an RCB star evolves from the 
merger of a He- WD and a CO- WD which has passed through a common envelope phase. He 
suggested that as the binary begins to coalesce because of the loss of angular momentum by 
gravitational wave radiation, the (lower mass) He- WD is disrupted. A fraction of the helium 
is accreted onto the CO- WD and starts to burn, while the remainder forms an extended 
envelope around the CO- WD. This structure, a star with a helium-burning shell sur rounded 
by a ~ 100 R Q hydrogen-deficient envelope, would resemble that of an RCB star ( IClayton 
19961 ). The merging times (~10 9 yr) might not be as long as previously thought, which makes 
the DD scenario an ap pealing altern ative to the FF scena rio for the formation of RCB stars 
rtlben et alJll997L Il996h . In addition. Isaio fc Jeffervi f|2002h suggested t hat a WD- WD merger 
could also account for the elemental abundances seen in RCB stars. iPandey et al.l (120061 ) 
have suggested a similar origin for the EHe stars. 

In order to judge the viability of this scenario in the light of our new observations, and 
in order to determine the direction of future more laborious and detailed modeling efforts we 
have investigated two separate questions in parallel. We have analyzed the thermodynamic 
conditions during the dynamic merger process of the two WDs (Section l5.ip . and we have 
investigated the necessary conditions to generate the observed abundance patterns from a 
nucleosynthetic point of view (Section [5. 2p . As we show in this section, these two approaches 
give compatible results, thereby boosting our confidence that the new observations will 
eventually be fully explainable in the WD merger scenario. 
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5.1. The Double Degenerate Merger 

The purpose of this section is to determine the thermodynamic conditions encountered 
during the dynamic merger process of a He- and a CO- WD. We discuss first the conditions 
that lead to complete, dynamical tidal disruption of the He-WD, which is the scenario 
in which the large overproduction of 18 may be possible. The critical parameter is the 
mass of the He-WD. After discussing the conditions for tidal disruption, we examine the 
thermodynamic conditions during disk accretion to determine if they allow the types of 
nucleosynthesis that might result in the observed overproduction of 18 O. 



5.1.1. Tidal Disruption 

To determine the extent and structure of the disk formed in the merger, we must examine 
the merger evolution. The merger process begins when the He-WD is brought so close to 
its companion (either by magnetic braking or gravitational wave emission) that it overfills 
its Roche radius. The orbit a l sepa ration at which a WD of this radius will overfill its Roche 



radius is given by lEggletonl (119831 ) : 



0.6g 2 / 3 + ln(l + q 1 / 3 ) 

A 0~RBe > C 1 ) 

where q = M He /M co with M co equal to the mass of the CO- WD and M He , Rn e are se t to the 
mass and radius of the He-WD, respectively. Assuming a T = 5/3 polytrope ap proximation 



of a W D equation of state, we derive an equation for the radius of the He-WD ( INauenberg 



1972|): 



R He = lO 4 (M He /O.7M )- 1 /3[i _ MHe/Mch] 1 / 2 ^^)" 5 / 3 ^, (2) 

where Mch ~ 1.4 is the Chandrasekhar mass and fi e is the mean molecular weight per 
electron of the WD. 

As matter flows from the He-WD onto a disk around its companion, both the Roche- 
overflow radius of the WD and the orbital separation of the binary change. The degenerate 
nature of the He-WD causes it to expand as it loses mass. If orbital angular momentum were 
conserved, the orbital separation would also increase as the low-mass He-WD accretes onto 
its heavier companion. The accretion process could then follow one of two possible paths: 
the He-WD radius increases faster than the orbital separation and the WD is disrupted over 
the course of one or two orbital periods, or the orbital separation increases faster than the 
WD radius and the He-WD accretes onto its companion on a timescale set by magnetic 
braking. To determine the path our objects follow, we compare the expansion of both the 
He-WD and the binary separation. 
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However, orbital angular momentum is not conserved; a fraction of it is transferred into 
angular momentum of the disk. Assu ming no mass is ejected, the angular mome ntum as a 
function of the evolving WD mass is jFrver et alill999l : [Podsiadlowski et al.lll992h : 

C l / 1\/T \ C2 

(3) 



A 



ML 



Mco 
M° co 



where G\ = — 2 + 2 jdisk and C2 = —2 — 27^1 and Jdisk is the sp ecific angular mom entum 
of the accretion disk material (see iFryer et al.lll999l . for details). iFryer et al.l (119991 ) found 
for the case of WD-black hole mergers that more than half of the angular momentum in the 
accreting material found its way into the disk or into the black hole spin (that is, jdisk ^ 0.5). 
It is likely the larger disks produced by these low-mass WD mergers will carry even more 
angular momentum. 

The orbital separation as a function of He-star mass for four separate binary systems 
is shown in Fig. [3J The solid curves show the path we expect the binaries to follow as the 
He-star loses mass (going right to left in decreasing He-star mass). The four curves denote 
four different values of jdisk- The dotted curve shows the separation for Roche- lobe overflow 
for this binary (again, a function of the He-star mass). There are two things to note in this 
figure. First, if jdisk > 0.5, the separation of the binary system will always be less than 
the separation for Roche- lobe overflow until the star is disrupted (~ 0.1 M Q ). This means 
that when the WD starts accreting, it will not stop until the companion WD is completely 
disrupted. The accretion process should proceed on a few orbit timescale: ~ 1000 s. Second, 
the orbital separation does not increase significantly during this accretion process for such 
high values of jdisk- Tidal effects will send so me material as far out as twice the initial 
separation (but very little will be ejected - see IFryer et al.l (Il999l )). But the bulk of the 



accreting material will not span more than twice the initial orbital separation at the onset 
of Roche-lobe overflow. 



5.1.2. Disk Accretion 



Based on WD merger calculations by lFryer et al.l (119991 ). we expect most of the matter 



in the He- WD to accrete onto its CO-WD companion. The He- WD is disrupted nearly 
instantaneously, forming a disk with a total radius of nearly 2 times the separation at the 
onset of Roche-lobe overflow. The accretion rate onto the CO-WD will definitely vary with 
time, but to obtain a preliminary estimate of it and of the temperature, we calculate the 
average value for these quantities. 

Under the a disk prescription for transporting angular momentum in the disk formed 
by the WD disruption, we can estimate the accretion time: t acc = Pdisk/« where -Pdisk is 
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the orbital period of the disk and a is the viscosity of the disk. Fig. H] shows this accretion 
time as a function of the mass of the CO- WD assuming a = 10~ 3 . The corresponding mass 
accretion rate (M acc ), roughly Mn e /t^ cc , is also plotted in Fig. HI The entire disk is likely to 
accrete in a few days, with accretion rates as high as 150M Q yr _1 . 

If this material accreted in a steady system, it would build up an atmosphere o n the WD. 
Accretion atmospheres have been studied i n detail for th e case of neutron stars (j Chevalier 
19891 ; iFryer et al.lll996l ; IColgate et al.lll993l ). iFryer et al.l (119961 ) derived the entropy, S, of 
this atmosphere as a function of the compact remnant mass, in this case, the WD mass 
(Mwd), position of shock (r s h oc k)E, and the accretion rate (M acc ). The pressure of this 
accreting material will be nearly dominated by photons. For a 7 = 4/3, radiation dominated 
gas, the entropy is given by: 



S = 27.9M^M-y 4 



p -3/8 
shock 



(4) 



where the Mwd is in M , M acc is in M & yr and r s h oc k is in 10 cm. As the atmosphere 
builds, the entropy at higher layers will decrease, and this negative entropy gradient will lead 
to rapid convection, producing a flatter entropy profile on timescales inversely proportional 
to the Brunt- Vaisala frequency, roughly 1-100 s. But the extent of this atmosphere, and 
hence, the lower limit on the entropy, is difficult to determine without actual calculations 
(indeed, given the dynamic nature of this accretion, such a simple picture may not even be 
correct). Typical initial entropies are roughly 10-30 erg K _1 for this accreting atmosphere, 
but it is possible that the convection can drive the value down by a factor of ~ 2. We can 
derive a temperature and a density of the atmosphere as a function of entropy: 



T = 2.262 x 10 8 5 ,-1 (i2/10 10 cin)- 1 K, 



(5) 



and 



p = 390(S/10k B nudeon -1 )- 4 (i2cra/10 9 cm)- 3 g cm" 3 . 



(6) 



The corresponding temperature for the atmosphere at the surface of the CO- WD lies in the 
1 — 2 x 10 8 K range if we assume the highest value for the entropy. The temperature could 
increase if mixing brings lower entropy material down, but is likely not to exceed 4xl0 8 K. 
The high temperatures and short dynamic nature of the merger process are supported by 
multi-dimensional SPH simulations as well (IGuerrero et al.l 120041 ). 



3 This is initially the CO- WD radius but it may move outward an order of magnitude 
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5.2. Implications from Nucleosynthesis 

Now that we have a rough physical picture of the merger process and the corresponding 
thermodynamic conditions we turn to the nucleosynthesis part of the problem. We ask what 
conditions are required for nucleosynthesis to produce the observed abundances. 



5.2.1. Observational Signatures ojB.dC and RCB Stars 



Typical RCB abundances are characterized by extreme hydrogen deficiency, en richment 
relatiye to Fe, of N, Al, N a, Si, S, Ni, the s-process elements and sometimes O (ILambert 
19861 ; lAsplund et al.l |2000| ) . Based on their abundance characteristics, RCB stars can be 



divided into a homogeneous majority and a diverse minority class, the latter distinguished 
by, among other things, lower metallicity, generally less strikin g H-deficiency, and extreme 
abundance ratios such as Si/Fe and S/Fe (lAsplund et al.ll2000l ). The isotopic carbon ratio, 
12 C/ 13 C, where measured, is very large (> 500) in most RCB stars. One exception i s the 
mino rity-class RCB, V CrA, where 12 C/ 13 C is estimated to be between 4 and 10 ((Rao 
20051 ). The elemental abundances of the HdC stars are sim ilar to those of the majority of 
RCB stars (jWarnerl 119671 ; lAsplund et al .11 1 9 9 71 : iKipperl 120021 ) . Lithiu m is seen in four of th e 
majority- class RCB stars and at least one HdC star, HD 148839 fiRao &: Lambertl Il996l ). 
Also, alt hough difficult to est imate, C/He is believed to be high in both classes of stars 
(-1-10%, lAsplund et aDboOot ). 



The surface compositions of HdC and RCB stars are extremely He-rich (mass fraction 
0.98), indicating that the surface material has been processed through H-burning. After 
H-burning via the CNO cycle, N is by far the most abundant of the CNO elements, because 
14 N has the smallest nuclear p-capture cross-section of any stable CNO isotope involved. 
However, the majority of RCB stars has \ogC/N = 0.30 and log N/O = 0.4 by number, 
equivalent to mass ratios of C/N = 1.7 and N/O = 2.2 jAsplund et al.lboooh . The N/O 
ratio represents the average for the majority RCB stars although the individual stars show a 
large scatter. Thus C is the most abundant and O the least abundant CNO element. These 
abundances are consistent if the material at the surfaces of HdC and RCB stars experienced 
a small amount of He-burning, as for example at the onset of a He-burning event that is 
quickly terminated. This partial He-burning would not significantly deplete He, but could 
be sufficient for some of the 14 N to be processed into 18 0. At the onset of He-burning, 13 C is 



4 There is still an unresolved C-problem in the spectral abundance analysis. The value, \ogC/N = 0.3, 
for the majority RCB stars is only true if the de rived C abundance r ather than the input C abundance is 
adopted, which rather would give log C/N = 0.9 (jAsplund et al.ll2000h . 
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the first a-capture reaction to be activated because of the large cross-section of 13 C(a, n) 16 0. 
Thus, a large 12 C/ 13 C ratio and enhanced s-process elements are both consistent with partial 
He-burning. 

5.2.2. Parametric Nucleosynthesis Model 

We constructed a parametric nucleosynthesis model for a quantitative estimate of the 
partial He-burning nucleosynthesis scenario (Fig.^J). We consider a single zone nucleosyn- 
thesis calculation with a solar initial abundance distribution. The density is assumed to be 
2 x 10 3 g cm -3 and the temperature initially is 5 x 10 7 K. We follow the nucleosynthesis with 
a nuclear network code until the H has been burned to a mass fraction < 10~ 4 . This situation 
is reached after approximately 125 yr. This first phase in our model represents H-burning, 
e.g., in the core or in a shell, that produces the He in the He- WD. The assumed density is 
unrealistically high, because, for simplicity we do not adjust the density in different phases 
of the model. In practice, varying the density will only affect the timescale associated with 
each burning phase. Here, the first phase serves only one purpose: to easily arrive at an 
abundance distribution that reflects that in a He- WD. 

At the end of the H-burning phase, we increase the temperature to switch smoothly 
into the He-burning regime (Fig.[3J. In our scenario, this temperature rise will occur as a 
result of the accretion of the He- WD onto the more massive CO-WD. Here we assume a 
temperature rise of 1% yr -1 . The timescale is not important here as everything scales with 
the density. The purpose is only to see how the relative abundances change as a function of 
temperature. 

We stop the increase when T = 1.65 x 10 8 K, which occurs at ~ 245 yr. For even slightly 
larger maximum temperatures, He-burning quickly reduces the abundances of He and 14 N to 
much smaller values than observed in the HdC and RCB stars. Lower temperatures require 
a much longer time for 14 N to be transformed into 18 O. From the merger considerations in 
the previous section, we anticipate that the period of He-burning cannot be very long if it 
is a result of the dynamical merging process. The exact timescale depends on the density 
evolution. After stopping the temperature increase, we allowed the network to continue for 
another 255 yr until abundances are reached that clearly disagree with those seen in the HdC 
and RCB stars. 



-14- 



5.2.3. Results 

Fig. [5] (upper right panel) shows the evolution of the dominant CNO isotopes. Initially 
the CNO H-burning equilibrium ratios, including the large 14 N abundance are established. 
Then, as the temperature increases, a modest increase in the 16 abundance can be seen 
(at 220 yr), resulting from the activation of 13 C(a,n) 16 0. Thereafter 16 remains largely 
unchanged because the temperature does not become high enough for 12 C(a, 7) 16 0. However, 
just after 13 C burning, 14 N begins undergoing a-capture, producing 18 0. This reaction very 
quickly makes 18 the most abundant of the CNO isotopes, a situation that lasts until 18 
begins capturing a particles to produce 22 Ne at which time it is overtaken by 12 C whose 
abundance has been increasing since t = 230 yr because of the triple-a reaction. 

The evolution of the elemental C/N and N/O ratios (lower left panel) is qualitatively 
consistent with nucleosynthesis after full completion of H-burning and at the beginning of He- 
burning. The N/O ratio decreases when the temperature is high enough to start He-burning, 
and eventually reaches the mean observed value at ~ 250 yr. The C/N ratio increases, but 
reaches the observed mean value somewhat later. However, both times are consistent with 
a very high and only slightly decreased He abundance at this early time of He-burning. 

The most pertinent information in Fig. 5 is the evolution of the 16 0/ 18 and 12 C/ 13 C 
ratios (lower right panel). The enormous production of 18 from 14 N leads to a rapid decrease 
of the 16 0/ 18 ratio at the onset of He-burning, which passes through the observed value 
for HD 137613 at around 25 yr. The 12 C/ 13 C ra tio increases to the point that it exceeds its 



observed lower limit of 500 (IFujita fc Tsujil 119771 ) at a slightly earlier time 



5.2.4- Admixture of H 

It is intriguing that most of the observed abundances of the HdC and RCB stars that 
have been discussed so far seem consistent with a nucleosynthetic origin at the onset of 
He-burning. One must keep in mind that this simulation consists of only a single zone, 
and that the temperature evolution is roughly assumed to demonstrate the various stages of 
nucleosynthesis in a DD merger. Nevertheless, one may complain that the two observations 
are not well matched. The C/N ratio increases too late, when 16 0/ 18 and N/O have already 
reached values far from those observed. The other inconsistency is the observed enhancement 
of s-process elements. In our model, the only source of neutrons is the small amount of 13 C 
from H-burning via the CNO cycle. This is insufficient by far to produce significant s-process 
enhancements. In addition, a good match for the CNO elemental and isotopic abundances is 
only obtained for a short period. Possibly the observed values are the result of a characteristic 
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termination timescale for this process and the associated nuc l eosynthesis. The adm ixture of 
H could lead to an enrichment of 14 N ( IRao fc Lambertlll996l ; lAsplund et al.ll2000l ). 



However, the CO- WD retains a small H-rich envelope, 1O _3 -1O _4 M depending on the 
WD mass. The accreting material will be shocked and heated, quickly leading to violent 
He-burning, which will induce convection. In addition, the accreted layer will be perturbed 
by the ongoing dynamic merging phase. All this will help to bring the original H-rich WD 
envelope up into the accreted He-burning material. In particular, the fact that the accreted 
elements have larger atomic masses make mixing of the H-rich CO-WD envelope into the 
accreted He-rich material inevitable. 

The question is when and how much H-rich material will be added. With our one-zone 
model we can only crudely explore the effect of such an H-admixture. In reality, there will 
be multi-zone mixing and burning at a range of densities and temperatures. Nevertheless, 
if the temperatures of the hottest layers are about right, one can obtain a zeroth-order 
approximation to what really happens. The admixture of a substantial fraction of the H-rich 
CO-WD envelope has to occur at a time when some 12 C has already been produced through 
He-burning, because a significant new source of 12 C is required for additional p-captures to 
produce the 13 C that can serve as a neutron source for the s-process. From Fig.[5l it is clear 
that triple-a reactions are only activated at a temperature of 1.6 x 10 s K at the assumed 
density. Then, the admixture must occur after that temperature has been reached. 

To test these ideas, we selected a CO-WD envelope mixture to be blended into the He- 
burning accreted material. The mixture comes from an actual post-AGB stellar evolution 
calculation with a core mass of 0.6 M & and an initial ZAMS mass of 2 M & . It has enhanced 
12 C (X c = 1.4 x 1CT 2 ) and 16 (X = 1.5 x 1CT 2 ), considerably larger than our initial choice 
of solar abundances (X c = 3 x 1CT 3 , X =8.3 x 1CT 3 ) due to several third dredge-up events 
during the TP-AGB evolution phase preceding the WD formation. The accreted He- WD 
material has a mass of order 0.1 M , while the CO-WD envelope has a mass of the order 
5 x 10 -4 Mq. Thus, we admix the envelope in a ratio 1:0.005 and renormalize to unity. We 
examine the effect of this admixture at the time when the 12 C abundance reaches a mass 
fraction of 5 x 10~ 3 . 

Fig. [6] shows the results of the network calculation. As the H-rich material is admixed at 
this early stage of He-burning, it encounters a large 18 abundance. Through two successive 
p-captures, 18 0(p, a) 15 N(p, a) 12 C, most of the 18 is transformed into 12 C. 19 F is formed 
by 18 0(p,7) 19 F and, because of the large abundance of 4 He also through 15 N(a,7) 19 F. 
Temporarily, the 19 F abundance exceeds the solar value by two orders of magnitude. We 
note that lines of F I were recently detected for the first time in visible spectra of a group 
of extreme He stars JPandev! lioOfl ). The 12 C nuclei, both freshly produced from triple-a 
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(primary), and coming from 18 make 13 C in turn producing neutrons and 16 through 
a-capture: 12 C(p, 7 ) 13 N(/3+) 13 C(a, n) 16 0. 

Due to the burning of 13 C, the 16 abundance is increasing quickly. The abundance of 
14 N also increases because of 13 C(p,7) 14 N. The admixture takes place just after the C/N 
ratio has increased significantly above the observed mean value. Due to the burning of 12 C 
and the creation of N, the C/N ratio drops below the observed value again. The opposite 
happens to the N/0 ratio. At the time of the admixture, 18 dominates the elemental 
abundance. Its quick decrease plus the increase in 14 N lead to a rapid increase in the N/O 
ratio above the observed ratio. 

The effects of the admixture on the 12 C/ 13 C and 16 0/ 18 ratios are dramatic. Due to 
the burning of 18 and the increase in 16 from the 13 C(a, n) 16 reaction, the 16 0/ 18 ratio 
quickly rises beyond the observed values. The 12 C/ 13 C ratio decreases, partly due to the 
admixture of 13 C from the CO-WD envelope (the first jump down to log 12 C/ 13 C « 4) and 
then partly due to the production of 13 C from 12 C(p,7). 

The effect on s-process elemental abundances should also be profound. Due to the 
production of the 13 C at high He-burning temperatures, neutrons are rapidly released. In 
this particular admixture test, we find r cxp = 0.13 mbarn -1 and a peak neutron density of 
N D > 10 11 cm~ 3 . T hus, the neutron expo sure is in the range that could explain the observed 



s-process elements ( lAsplund et al.ll2000i ). However, the neutron density is far too high to 



produce an s-process signature. 

In summary, it seems that the admixture of H can aid in solving the delay in the match of 
the elemental and isotopic ratios. It also has, in general, the effect of reversing the effect of too 
much He-burning. For example, the He-burning leads to the decrease of 16 0/ 18 to values far 
below the lowest observed. H-admixture has the effect of reversing this trend and increasing 
the O isotopic ratio again. The same can be observed for the other abundance ratios. We 
therefore speculate that H-admixture in a realistic, highly mixed accretion environment can 
help to stabilize the nucleosynthesis results leading to the significant coherence of abundance 
patterns observed in HdC and RCB stars. 



6. Conclusions 

Observations of the first overtone CO bands in cool HdC and RCB stars have resulted 
in two significant findings: (1) in all HdC stars and most or all RCB stars where oxygen 
isotopic ratos can be measured, 18 O is enhanced relative to 16 O by factors of 100-1000 
compared to standard Galactic values; and (2) the lowest 16 0/ 18 ratios are found in the 
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HdC stars. Our models show that nucleosynthesis in the dynamically accreting material of 
CO-WD-He-WD mergers, previously suggested as the source of RCB stars, may provide 
a suitable environment for significant production of 18 and the very low 16 0/ 18 values 
as observed. The suggestion that FF events produce RCB stars appears to be weakened 
due to a mismatch of oxygen (this paper) and carbon isotopic abundances, as well as by 
the apparent inconsistency in timescale of the RCB-mimi cking portion of the FF episode 
and that required to account for the number of RCB stars (jClayton et al.ll2006l ). Therefore, 
we have explored the idea that HdC and RCB stars originate in the mergers of CO- and 
He-WDs. 

In the WD merger, it is partial He-burning at just the right temperature and duration 
that is the likely generator of small 16 0/ 18 ratios. In addition, the high values of 12 C/ 13 C 
observed in HdC and RCB stars can probably be reproduced at the same time. The observed 
s-process abundances in these stars as well as the delayed increase of the C/N ratio to the 
observed mean value motivated the investigation of the effects of admixing H-rich material. 
This material would, in the DD merger /accretion scenario, be the thin CO-WD envelope. 
We found that qualitatively the instantaneous admixture of H-rich material resets all ob- 
servables back toward the observed values. The fact that all of the H-rich material is added 
instantaneously leads to artefacts such as the extremely large neutron density, as well as the 
very small 12 C/ 13 C ratio. In reality, the admixture will take place over a finite time and in 
smaller quantities. A slower rate of admixing could produce the observed abundances. 

It is interesting to no te that some RCB stars show enhan ced Li abund ances, as does the 
FF sta r Sakurai's Object jLambertlll986l : lAsplund et al. Il998h . As shown by lHerwig fc Langerj 



(120011 ). Li enhancements are consistent with the FF scenario. However, the production of 
18 requires temperatures large enough to at least marginally activate the 14 N(a, 7) reaction. 
The a capture on Li is eight orders of magnitude more effective than on 14 N. For that reason 
the simultaneous enrichment of Li and 18 is not expected in the WD merger scenario. This 
is an important argument against the FF evolution scenario as a progenitor of RCB and 
HdC stars with excess of 18 0. The presence of Li in some RCB stars, as well as in Sakurai's 
Object, is extremely difficult to explain with a WD merger. Since 18 strongly supports 
the DD merger /accretion scenario for most of the stars studied so far, the obvious solution 
is that there are (at least) two evolutionary channels leading to RCB, HdC and EHe stars, 
perhaps with the DD being the dominant mechanism. Unfortunately the division between 
majority- and minority- class RCB stars does not lend itself naturally to t his explanation, 
since Li has only been detected in the majority group (lAsplund et al.ll2000l ). 



Contrary to the work of lSaio fc Jefferyl (120021 ) , who assumed an accretion rate of 10 5 M / yr 
to model WD mergers with a hydrostatic evolution code, our estimate of the merger process 
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leads to a picture of a rapid, dynamic event that only lasts a few days, and during which 
high temperatures can be obtained. In the merger scenario, rapidity appears to be essential 
to producing the conditions required to account for our observations of extremely abundant 
18 0. 

Understanding the RCB and HdC stars is a key test for any theory that aims to explain 
hydrogen deficiency in post-AGB stars. Confirmation of the WD merger scenario will allow 
the use of RCB and HdC stars as probes for WD merger simulations. The ability to model 
the rates of these low-mass WD mergers will help us to underst and the rates of more massive 
mergers that make type la supernovae (IBelczynski et al.ll2005l ). 



Most of these observations were obtained during program GS-2005B-Q-20 at the Gemini 
Observatory, which is operated by AURA under a cooperative agreement with the NSF on 
behalf of the Gemini partnership of Argentina, Australia, Brazil, Canada, Chile, the United 
Kingdom and the United States of America. We are indebted to the Gemini staff for its 
expertise in obtaining these data. We also thank the staff of UKIRT, which is operated by 
the Joint Astronomy Centre on behalf of the U.K. Particle Physics and Astronomy Research 
Council. TRG's research is supported by the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., on behalf of the international 
Gemini Partnership. This work was carried out in part under the auspices of the National 
Nuclear Security Administration of the U.S. Department of Energy at Los Alamos National 
Laboratory under Contract No. DE-AC52-06NA25396, and funded by the LDRD program 
(20060357ER). This work is LANL report LA-UR-07-0198. This work was carried out in 
part in collaboration with the NSF funded Joint Institute for Nuclear Astrophysics (JINA). 
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Table 1. CO Bandhead Absorption Depths in RCB and HdC Stars' 



Star 


Tel. 


Exp. 




12 G 


16Q 




12 C 


18Q 








Time (s) 


2- 





3-1 


2- 





3- 


1 


HD 175893 


GS 


480 





06 


0.07 





15 





11 


HD 182040 


GS 


56 





03 


0.05 





05 





06 


HD 137613 


UKIRT 


480 





12 


0.20 





18 





21 


WX CrA 


GS 


480 





18 


0.20 





15 





13 


S Aps 


GS 


160 





54 


0.51 





26 





15 


SV Sge 


GS 


80 





51 


0.47 





24 





11 


ES Aql 


GS 


480 





35 


0.31 





09 





07 


Z Umi 6 


Bok 


60 





34 


0.34 





07 





02 


U Aqr 


GS 


720 





16 


0.24 





03 


<0 


02 


HD 148839 


GS 


480 


<0 


05 


<0.03 


<0 


03 


<0 


02 


HD 173409 


GS 


960 


<0 


03 


<0.02 


<0 


02 


<0 


02 



a Based on flux ratios at rest vacuum wavelengths 2.2930 and 
2.2940 fim, 2.3215 and 2.3235 //m, 2.3480 and 2.3495 /im, and 2.3770 
and 2.3790 /im for 12 C 16 2-0 and 3-1 bands and 12 C 18 2-0 and 
3-1 bands, respectively. 



b Spectrum originally published by lTenenbaum et al.l (120051 ). 
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Table 2. Estimated Isotopic Oxygen Abundances 



Star 


HdC/RCB 


T e// 


abs(C 16 0)/abs(C 18 0) 


16 0/ 18 a 


HD 175893 


HdC 


5500 


0.4 


0.2 


HD 182040 


HdC 


5600 


0.6 


0.3 


HD137613 


HdC 


5400 


0.7 


0.5 


WX CrA 


RCB 


5300 


1.3 


1 


S Aps 


RCB 


5400 


2.5 


4 


SV Sge 


RCB 


4000 


2.5 


4 


ES Aql 


RCB 


5000 


3 


6 


Z UMi 


RCB 


5000 


>4 


>8 


U Aqr 


RCB 


6000 


>6 


>12 


HD 148839 


HdC 


6500 






HD 173409 


HdC 


6100 







Estimated uncertainty is a factor of two. 



References. — - lAsplund et al.l (119971 ); 
Bergeat et all (boOlh 



Lawson et al, 



(11990( 1: 
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2.25 2.3 2.35 2.4 2.45 2.5 

Wavelength (/urn) 



Fig. 1.— UKIRT spectrum of HD 137613. The wavelengths of the bandheads of 12 C 16 and 
12 C 18 are indicated. 
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Fig. 2.— 2.28-2.45/im spectra of RCB and HdC stars, with wavelengths of 12 C 16 and 
i2q18q b anc Lheads indicated by vertical lines. The stars are ordered by the strengths of the 
i2q18q k anc i s _ g rs t three spectra at the top and the last two at the bottom are of HdC 
stars. The rest are of RCB stars. 
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Fig. 3. — Evolution of a WD binary as a function of the masses of the two stellar compo- 
nents. The start of the evolution is where all the lines meet at a point (right side of the 
plots). The solid lines show the evolution of the orbital separation of the binary as mass 
transfers from the He- WD onto the CO- WD. The different curves correspond to varying 
amounts of angular momentum in the disk as a function of the mean binary angular mo- 
mentum: 0.2 means very little angular momentum goes into the disk, most stays in the 
binary whereas 0.8 means that the accreting material retains most of its angula r momentum 
as it for ms an accretion di sk. Models of WD/BH mergers (IFryer et al.lll999l ) and NS/BH 
mergers (IJanka et al.lll999l ) all suggest that this number is at least ~ 0.5. The dotted line 
denotes the binary separation below which Roche-lobe overflow occurs. Because degeneracy 
dominates the pressure of the He- WD, as it loses mass, it expands, and the separation at 
which Roche-lobe overflow occurs also increases. We start the evolution with the orbital 
separation set to the Roche overflow separation (the onset of Roche- lobe overflow). If the 
subsequent evolution causes the separation to drop below the Roche overflow separation, 
the accretion will develop dynamically and will lead to the disruption of the He- WD on an 
orbital timescale. The four panels correspond to two CO- WD masses (top - 0.9 M , bottom 
- 0.6 M ) and two He- WD masses (left - 0.4 M , right - 0.2 M ). 
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Fig. 4. — Properties of the accretion atmosphere as a function of the C/O mass. The solid and 
dotted curves correspond to 0.4 and 0.2 M Q He stars, respectively. The top left panel shows 
the accretion timescale, the bottom left panel shows the corresponding accretion rate. The 
largest uncertainty for both of these terms is the choice for the viscosity (we have assumed 
a pa 10 -3 ). This value could be a factor of 100 larger, producing accretion timescales a factor 
of 100 shorter and accretion rates a factor of 100 larger. If we assume this accretion builds 
an equilibrium atmosphere on top of the CO- WD, we can then estimate the temperature 
and density at the base of this atmosphere. This calculation depends upon the size of the 
atmosphere and also the accretion rate. The latter is the biggest uncertainty and could lead 
to an increase in the temperature of, at most, a factor of 2, and an increase in the density 
of roughly a factor of 10. We note, however, that nucleosynthetic constraints argue for our 
standard calculation. 
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Fig. 5. — Time evolution. Upper left: H and He as well as temperature. Upper right: CNO 
isotopes, and temperature. Lower left: C/N and N/O ratios and mean observed ratios of 
majority of RCB stars. Lower right: O and C isotopic ratios and observed ratios (straight 
lines) . 
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Fig. 6. — Upper left: Same as Fig. [5] for the short period around the time of admixing the 
H-rich CO- WD envelope material to the He-burning accreting material. Time evolution of 
CNO isotopes. Each mark along the 13 C line represents a time step in the network evolution. 
Upper right: Time evolution C/N and N/O ratios and mean observed ratios of majority of 
RCB stars. Lower left: Time evolution of O and C isotopic ratios and observed ratios 
(straight lines). 



